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STABLE  CRACK  GROWTH  IN  ALUMINUM  TENSILE  SPECIMENS 

A  AAA 

B.S.-J.  Kang  ,  A.S.  Kobayashi  and  D.  Post 

ABSTRACT 

Post's  white  light  moire  interferometry  was  used  to  obtain  sequential 
records  of  the  transient  u^-displacement  fields  associated  with  stable  crack 
growth  in  7075-T6  and  2024-0,  single  edge  notched  (SEN)  specimens  with  fatigued 
cracks.  The  u^-displacement  fields  were  used  to  evaluate  the  crack  tip  opening 
displacement  (CTOD),  far  and  near-field  J-integral  values,  Dugdale  strip  yield 
model,  William's  polynomial  function  and  the  HRR  fields. 

INTRODUCTION 

Crack  growth  in  ductile  material  can  be  divided  into  three  stages,  namely, 
1)  plastic  yielding  and  the  onset  of  stable  crack  growth,  2)  stable  crack 
growth  and  3)  rapid  tearing.  Since  the  measured  crack  velocity  during  rapid 
tearing  is  less  than  5  percent  of  the  dilatational  wave  velocity  Cl, 2, 31,  rapid 
tearing  and  stable  crack  growth  can  be  considered  as  quasi-static  deformation 
processes.  The  crack-tip  state  for  rapidly  tearing  and  stably  growing  cracks, 
however,  are  different  from  that  of  a  stationary  crack.  Asymptotic  analyses  of 
a  stationary  crack  in  an  elastic,  perfectly  plastic  solid  under  infinitesimal 
deformation  3how  that  the  strains  vary  as  1/r  but  for  a  growing  crack  the 
strains  vary  as  ln(l/r)  C4-81.  On  the  other  hand,  numerical  studies  C9-131  on 
stable  crack  growth  do  not  address  the  crack  tip  singularity  problem  but 
discuss  the  somewhat  near  stress  field  surrounding  the  crack  tip. 

Numerous  fracture  parameters  which  characterize  stable  crack  growth  under 
small-scale  yielding  condition,  such  as  average  crack  opening  angle  (COA)  C141, 
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crack  tip  opening  angle  (CTOA)  1153,  crack  tip  opening  displacement  (CTOD) 
1161,  critical  strain  117,181,  energy  release  rate  1191,  crack  tip  force  1203, 
J-resistance  curve  C211  and  tearing  modulus  (T)  1221,  have  been  proposed.  Of 
these,  the  crack  tip  opening  angle  (CTOA)  or  displacement  (CTOD)  was  3hown  to 
be  suited  for  modeling  stable  crack  growth  and  instability  during  the  fracture 
process  115,20,211. 

Under  large-scale  yielding  condition,  however,  there  is  no  analytical 
solution  available  for  stable  crack  growth.  Attempts  have  been  made  to  extend 
the  ductile  fracture  criteria  for  small-scale  yielding  and  stable  crack  growth 
to  large-scale  yielding.  The  few  results  published  to  date  11,2,15,201 
indicate  that  under  limited  conditions,  the  CTOA  or  CTOD,  are  plausible  ductile 
fracture  criteria.  The  purpose  of  this  paper  is  to  present  preliminary 
experimental  findings  on  the  crack  tip  parameters  which  control  the  initiation 
and  propagation  of  stable  crack  growth.  An  approximate  J-integral  evaluation 
procedure  based  on  u^-displacement  field  is  also  presented. 

ANALYTICAL  BACKGROUND 
( 1 )  J-integral 

For  two-dimensional  problems  of  materials  governed  by  nonlinear  elasticity 
and  deformation  plasticity  theory  subjected  to  monotonically  loading  condition, 
the  J-integral  is  defined  as  1231 

J  =  W  dy  -  ?  •  §£  ds  (1) 

p  oX 

where 

F  :  contour  surrounding  the  crack  tip 
f  :  traction  vector  along  the  contour 
u  :  displacement  vector  on  the  contour 
W  :  strain  energy  density  on  the  contour 


In  the  following,  an  experimental  procedure  for  direct  far-field  and 
approximate  near-field  J-integral  measurement  is  introduced.  The  underlining 


methodology  is  to  compute  the  J  value  with  only  the  u^-displacement  field  which 
is  obtained  from  a  single  moire  interferometry  recording. 


Far-field  J-inteqral  Measurement 
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where  a0  is  the  yield  stress,  N  is  the  strain  hardening  constant  and  a  is  a 
dimensionless  material  constant,  a  is  calculated  from  Equation  (2c)  tor  a 

yy 

given  e  which  can  be  determined  experimentally  from  the  moire  data. 

Thus  the  integral  value  of  Equation  (1)  along  the  vertical  edges  of 
segments  12  and  34»  is 
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where  i  is  the  ith  segment  of  the  contour. 

Along  the  horizontal  segment  23,  dy  is  zero  and  the  first  term  of  the 

integrand  in  Equation  (1)  vanishes.  The  traction,  ?,  along  this  segment  are 

T  =  a  and  T  =  t  .  At  this  point  we  assume  that  the  shear  stress,  t  , 
y  yy  x  xy  r  xy 

and  the  displacement,  ux  ,  are  negligible  along  segment  23,  This  assumption 
is  justified  if  segment  23  is  sufficiently  far  away  from  the  crack.  The 
integral  value  of  Equation  (1)  along  segment  23  thus  becomes 
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Aqain,  for  the  a  term,  the  same  stress-strain  relation,  e.g..  Equations  (2) 
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is  used.  Finally,  the  J-integral  value  is  given  by 


J  —  2(J  +  J,  ) 


(5) 


The  above  experimental  procedure  for  determining  J-integral  value  was 
carried  out  using  strain  gayes  and  linear  variable  displacement  transducers  at 
discrete  points  along  the  specimen  boundary  C24-261.  Since  the  test  data  in 
these  references  were  obtained  from  few  locations.  Equation  (5)  could  only  be 
evaluated  at  a  few  discrete  locations.  Moire  interferometry,  on  the  other 
hand,  provides  an  easy  alternative  for  implementing  this  method  with  better 
accuracy.  Since  it  yields  highly  sensitive  displacement  field,  the  approximate 
analysis  proposed  above  requires  only  a  single  u^-displacement  moire  field  for 
calculating  the  J-integral. 


Near-tield  J-integral  Measurement 

Whil-  the  above  procedure  is  valid  tor  far-field  J-integral  evaluation, 
its  validity  for  the  near  field  integration  contour,  such  as  the  inside 
rectangular  contour  shown  in  Figure  1,  must  be  justified.  First,  we  will  show 
that  the  above  far  field  J-integral  measurement  procedure  is  a  reasonable 
approximation  for  the  near-field  J  value  in  a  linear  elastic  field. 

Consider  a  rectangular  contour  around  the  crack  tip  as  shown  in  the  legend 
of  Figure  2.  For  a  linearly  elastic  material,  the  J  integral  along  the 
horizontal  segment  13  can  be  expressed  in  terms  of  displacements  u^  and  u^  as 
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where  G  is  the  shear  modulus,  v  is  the  Poisson's  ratio  and 
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The  second  term  of  the  integrand  in  Equations  (ba)  and  (fob)  were  neglected 
in  evaluating  the  far-tield  J-integral.  For  the  near-field  J-integral 
evaluation,  we  will  also  assume  that  the  contour  integrals,  as  represented  by 
Equation  (3)  and  (4),  can  be  used.  The  error  due  to  such  assumption  is 
evaluated  in  the  following. 

Consider  a  crack,  in  a  two  dimensional  linear  elastic  material.  The  mode  I 
crack  -tip  displacements  are 
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where  r  and  6  are  the  polar  coordinates  with  the  origin  at  the  :rack-tip, 
is  the  mode  1  stress  intensity  factor  and 
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\  plane-  strain . 
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After  substituting  Equations  (7a)  and  (7b)  into  Equations  (6a)  and  (6b),  the 
first  and  second  terms  of  the  integrand  in  Equations  (6a)  and  (6b)  are 
evaluated  along  a  non-dimensionalized  half  square  contour,  01234 ,  as  shown  in 


r 
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Figure  2.  Results  of  the  numerical  integration  using  either  the  first  and 
second  terms  or  the  first  term  alone  in  Equations  (6a)  and  (6b)  as  one 
traverses  along  the  half  contour  are  plotted  in  Figures  2  and  3.  In  these 
figures,  the  former  and  latter  J-integral  values  are  denoted  as  "theoretical" 
and  "approximate"  values  respectively.  Notable  is  the  close  proximity  between 
the  theoretical  and  approximate  summation  of  AJ,  £AJ,  along  the  contour  before 
entering  segment  34.  The  nondimensionalized  J=EAJ  values  at  point  4  shows 
about  14%  difference  between  the  theoretical  and  approximate  IAJ  values. 

We  further  evaluate  the  validity  of  this  J  approximation  procedure  tor  a 
crack  tip  region  characterized  by  a  Hutchinson-Rice-Rosengren  singular  field 
[27,281.  Again  consider  a  rectangular  contour  surrounding  a  crack  tip  as  shown 
in  the  legend  of  Figure  4.  The  HRR  stress,  strain  and  displacement  field 
within  this  rectangular  region  can  be  expressed  as  [27,281 
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where  I  is  a  dimensionless 
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of  0.  For  the  approximate  and  as  represented  by  Equations  (3)  and  (4), 
the  needed  a  and  W  can  be  represented  as 
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Equation  (9a)  and  (9b)  represents  the  plastic  components  of  Equations  (2c)  and 
(2d)  where  the  elastic  components  are  assumed  negligible  in  the  region 
characterized  by  the  HRR  field. 

The  approximate  J  can  then  be  evaluated  by  substituting  Equations  (8a)  and 
\8b)  into  Equation  (9)  and  evaluating  the  integral  along  the 
non-diraensionalized  half  square  contour.  Also,  the  theoretical  J  is  evaluated 
by  substituting  Equation  (8)  into  Equation  (1)  and  evaluating  the  integral 
along  the  non-dimensionalized  half  square  contour.  A  state  of  plane  stress 
with  N-2,  5,  50  were  chosen  for  this  analysis.  Numerical  values  for 
%  \ 

a.  .(0),  e.  (0)  and  u.(0)  were  obtained  from  [291.  The  results  are  shown  in 
1J  ij  i 

Figures  4,  5  and  6.  Good  agreement  between  the  theoretical  and  approximate 
J- integral  are  noted. 

The  results  of  Figures  2  through  6  suggest  that  the  approximate  J  as 
determined  by  the  far-field  solution,  is  reasonably  correct  when  used  in  a  HRR 
dominated  crack  tip  region.  However,  when  used  in  a  crack  tip  field  dominated 
by  linear  elasticity,  the  error  is  noticeable.  Figures  2  and  3  show  that  this 
error  is  generated  during  the  last  integration  path  or  along  the  vertical 
contour,  line  34,  indicating  that  the  assumed  uniaxial  tension  state  is  not  a 
reasonatle  approxima t ion  ot  the  true  state  of  elastic  stresses  along  line  34. 
in  contrast,  both  the  assumed  uniaxial  tension  state  ot  stress  and  the  t:ue 
state  ot  stresses  ot  the  HER  field  along  line  34  .have  negligible  ettoct  or.  the 


J  value  as  evidenced  in  Figures  4,  5  and  6  with  the  flat  portion  along  the  last 
integration  path  or  line  34,  and  thus  the  approximate  procedure  of  evaluating  J 
works  reasonably  well.  This  induced  error  in  the  elastic  crack  tip  stress 
field  can  be  reduced  it  line  34  is  situated  within  the  region  of  uniaxial 
tension  or  more  specifically  along  a  tree  boundary.  As  will  be  shown  later, 
under  such  restriction  the  approximate  J  will  provide  reasonably  accurate  J 
values  in  an  elastic  fracture  specimen. 

The  J  value  can  also  be  linked  to  the  crack  tip  opening  displacement 
(CTOD)  through  C30J  as 

6t  =  J/a0  (10) 

where  Cl,  values  can  be  found  in  [293. 


(2)Duqdale-Barenblatt  Strip  Yield  Model 

For  an  elastic  perfect-plastic  mp’-erial,  the  crack-tip  displacement 
fields  tor  a  Mode  I  plane  stress  IXigdale-Barenblatt  strip  yield  model  can  be 
expressed  C313  as 
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and  K  ^  is  the  Mode  I  stress  intensity  factor. 

The  CTOD  tor  the  Dugdale-Barenbla tt  strip  yield  model  become  [31 J 


Or 


1  ° 

CTOD  =  :~(4-4v)  —  r  , 
2G  v  y\ 


r  -  r 
X 


r  +  r 


-(-* - )  log  (-* - )V 

r  2r  3  r  -  r  ) 

'  i/  v  vr 


+  “  I  j(-l)n+1d  rn  1/ZF  ( n,Tt,  v) -M -1  )d  r  F  (n,ir,v)i 
2b  n-il  zn-i  3  2nni  ) 


( 12) 


EXPERIMENTAL  APPROACH 

White  light  moire  interferometry  L32J  was  used  to  obtain  a  single-trame 
record  of  static  and  dynamic  displacement  fields  surrounding  the  crack  tip  in 
slowly  and  rapidly  fracturing  7u7F-Tb  and  2U24-0  aluminum  2EN  specimens. 
Figures  7  and  d  show  the  optical  system  which  utilizes  a  compensator  grating  ot 


halt  frequency,  f/2,  where  f 


1200  lines/mra,  to  illuminate  the  reference  and 


specimen  gratings  ot  tuil  and  hair  frequencies,  respectively.  The  achromatic 
light  emerges  from  the  compensator  as  monochromatic  xignt  beams  at  ditterent 
diffraction  angles  and  generates  the  same  moire  pattern  tor  each  wave  length. 
The  camera  records  the  scalar  sum  ot  the  light  intensities  associated  with 
various  wave  length  and  thus  much  or  the  original  white  light  intensity  is 
recovered.  When  an  incoherent  light  source  is  used,  the  gap  between  the 
reference  and  active  gratings  must  be  small.  This  white  light  moire 
ir.u-r  tercmet :  y  provides  the  high  sensitivity  associated  with  high  frequency 
o  c  iru.  >:.i  bright  light  source  using  a  relatively  simple  experimental 

.  i  gn*  m^itv  rnnge  patterns  were  recorded  on  a  35  mm  camera 
w.  •  ■  i  ’  :  -j!  filtering.  A  motor-driven  camera  provided  up  to  b 
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loadings.  Tne  specimen  c  r.t  i  rvraM  : r.,  ;:.,c  -1  .  properties  and  the  two  material 
coefficients  f  or  the  power  Hardening  sir  •  •ss- strain  1  eiati-nx  are  shown  in 
Figure  These  material  pr  , per  ties  indicate  tnat  aluminum  'u/5  Tb  is 

essentially  an  elastic  -  per  let;*,  p.a.ti:  materia  j  while  ..014  •  is  1  strain 

hardening  material.  Figures  iu  and  11  ire  typical  wh ;  f  e  light  moire 
interferometry  fringe  ra » t  erne  ■■1  aluminum  '"/h  T»_  and  30x4-0  HEN  specimen  with 
star;-'  Vick  u  wr..  Th<-  if  pr  ./  ••  ••v.i.n.i'  i  n  pt  ocedur  e  was  used  to 
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Since  7o75-Tt>  aluminum  is  a  relatively  brittle  material  with  the  crack-tip 
being  surrounded  with  small  scale  yielding,  the  tar-tield  J  value  as  well  as 
the  near-rield  J  value  outside  ot  the  yield  zone  can  be  determined  by  elastic 
analysis.  These  elastic  values  are  used  to  verity  the  accuracy  or  the 
experimental  and  data  reduction  procedures  used  in  this  paper. 

Figure  12  shows  the  log-log  plot  ot  the  u^-displacement  versus  radial 
distance  up  to  marked  boundary  where  the  slope  was  0.5+0.05.  Experimental 
deviation  m  the  slope  of  log  u^  versus  log  r  curve  was  determined  by  linear 
regression  ot  a  straight  line  fitting  through  the  data  points  in  Figure  12  and 
then  computing  the  percentage  deviation  in  slope  from  the  crack  tip.  The 
average  slope  of  1/2  in  the  vicinity  of  the  crack  tip  indicates  that  the 
elastic  field  prevailed  in  this  specimen  up  to  a  distance  1.2  mm  from  the 
crack  tip.  The  approximate  J  values  which  were  determined  by  the  above 
mentioned  J  evaluation  procedure  are  shown  in  Table  1.  Also  shown  for 
comparison  purpose  in  Table  1  are  the  corresponding  stress  intensity  factor 
values  computed  by  K  =  JJ*E  using  the  J  values  obtained  from  the  moire 
fringes  and  the  K  values  computed  by  using  the  formula  in  ASTM  STP  410  C33J, 
the  William  s  polynomial  function  C343  and  the  Dugdale  strip  yield  model. 
Figure  13  is  the  corresponding  plots  of  the  stress  intensity  factor,  K,  versus 
applied  load.  Good  agreements  between  the  measured  K  and  that  computed  by  ASTM 
STP  410  results  are  noted.  Also  shown  in  Table  1  are  the  experimentally 
measured  and  the  computed  CTOD  values  based  on  the  Dugdale-Barenblat t  strip 
yield  model.  The  computed  CTOD  value  were  obtained  by  least  square  fitting 
Equation  (11)  with  n=2  which  is  a  four  parameter  characterization  of  the  crack 
tip  stress  field  to  the  u  -displacement  field  of  the  moire  tnnqes.  The 

y 


parameters,  (e.g.  d  and  r  )  were  then  back  substituted  to  Equation  (12)  and 

n  y 

the  CTOD  value  was  computed.  Figure  14  shows  CTOD  plots  ot  the  7075-T6  fatigue 
precracked  specimens  versus  applied  load.  The  same  CTOD  values  were  observed 
before  the  onset  of  unstable  crack  growth  in  this  fatigue  pre-cracked  SEN 
specimens. 

Table  2  shows  the  approximate  J  values  which  were  determined  along  the 
three  contours  in  the  7075-Tb  SEN  specimen  shown  in  Figure  10,  for  ten 
sequential  moire  fringe  patterns  ot  stable  crack  growth.  As  expected  in  this 
elastic  specimen,  the  J-values  along  these  three  contours,  tar  to  near  field 
contours  are  in  good  agreement  with  each  other. 

2024-0  SEN  Specimens 

A  2024-j  SEN  specimens,  with  fatigue  precrack,  were  tested  to  failure. 
Unlike  the  7075-Tb  specimens,  no  unstable  crack  growth  were  observed  in  these 
specimens  which  exhibited  large  scale  yielding.  Figure  15  shows  the  log-log 
plots  ot  the  Uy-displacement  versus  distance  to  the  crack  tip  ot  Figure  11. 
The  average  slopes  near  the  crack  tip  within  the  marked  region  in  Figure  15 
is  l/b_f  0.U2  which  is  the  predicted  exponent  for  a  HRR  displacement  fields 
L27,2bJ.  Table  3  shows  the  approximate  J  values  obtained  from  the  sequential 
moire  interferometry  recordings  or  the  test  specimen  along  three  different 
paths  as  shown  in  Figure  11  tor  each  frame.  These  results  show  that  within 
the  relatively  short  crack  extension  ot  0.75  mm,  J  is  still  a  valid  parameter 
tor  characterizing  the  crack  tip  L35J  and  Equation  (10)  is  valid  within 
this  loading  range.  The  path  independency  ot  the  measured  J  values  tor  each 
I rame  is  an  experimental  validation  ot  the  J  estimation  procedure  proposed  in 
this  paper.  Figure  lb  shows  the  increases  in  the  approximate  J  values,  which 


are 


consistent  with  published  results  120,211,  with  crack  extension  tor 


the  2024-0  specimen.  Notable  m  Figure  16  is  the  slopes  dJ/da  of  the  J 
values  which  remain  constant  during  the  initial  short  crack  extension  ot 
0.4-0. 5  mm  and  its  continuous  decrease  beyond  this  short  crack  extension. 
Figure  17  shows  the  increases  in  CTOD  values  with  crack  extension.  These  CTOD 
values  can  be  linked  to  the  corresponding  measured  J  values  through  Equation 
1 10)  where  Figure  18  show  the  variations  ot  the  experimental  determined  D^ 
values  which  agree  favorably  with  the  value  determined  by  the  HRR  field, 

le.g.,  Dn=0.33  for  2024-0  aluminum  material). 

CONCLUSIONS 

1.  White  light  moire  technique  was  used  to  determine  the  u^-displacement 
field  of  stably  growing  cracks  in  7075-Tb  and  2024-0  aluminum  SEN 
specimens . 

2.  An  procedure  for  estimating  J  from  the  recorded  u^ -displacement  field  was 
developed.  This  approximate  J  agrees  reasonably  well  in  a  HRR  crack  tip 
field  but  required  special  handling  when  used  in  an  elastic  crack  tip 
field . 

3.  The  approximate  J  and  CTOD  at  the  onset  and  during  stable  crack  growth 
were  recorded.  Limited  data  suggests  that  there  exist  a  constant  CTOD  tor 
unstable  crack  propagation  m  7075-Tb  aluminum  SEN  specimen.  HRR  field 
dominates  stable  crack  growth  in  2024-0  aluminum  SEN  specimen. 
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Table  1  Measured  and  Calculated  J,  1C  and  CTOD  of  7075-T6  Aluminum  SEN 
Specimen  with  Fatigue  Precrack. 
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1*  :  based  on  ASTM  STP  410  K  evaluation  procedure. 

2*  :  based  on  J  evaluation  procedure,  e.g.  K=iJ*E. 

3*  :  based  on  the  William  s  polynomial  function. 

4*  :  based  on  the  Dugdale-Barenblatt  strip  yield  model. 

**  :  based  on  the  Dugdale-Barenblatt  strip  yield  model. 

***:  rapid  crack  growth. 


Table  2  Measured  Approximate  J  Values  for  Different  Contours  in 
7075-T6  Aluminum  SEN  Specimen  with  Fatigued  Precrack. 
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Table  3  Measured  Approximate  J  Values  for  Different  Contours  in 
2024-0  Aluminum  SEN  Specimen  with  Fatigued  Precrack. 
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